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ABSTRACT 

We explore the possibility that the steepening observed shortward of 1000 A 
in the energy distribution of quasars may result from absorption by dust, being 
either intrinsic to the quasar environment or intergalactic. We find that a dust 
extinction curve consisting of nanodiamonds, composed of terrestrial cubic dia- 
monds or with surface impurities as found in carbonaceous chondrite meteorites, 
such as Allende, is successful in reproducing the sharp break observed. The inter- 
galactic dust model is partially successful in explaining the shape of the composite 
energy distribution, but must be discarded in the end, as the amount of crys- 
talline dust required is unreasonable and would imply an improbable fine tuning 
among the dust formation processes. The alternative intrinsic dust model re- 
quires a mixture of both cubic diamonds and Allende nanodiamonds and provide 
a better fit of the UV break. The gas column densities implied are of the order 
10 20 cm -2 , assuming solar metallicity for carbon and full depletion of carbon into 
dust. The absorption only occurs in the ultraviolet and is totally negligible in the 
visible. The minimum dust mass required is of the order ~ O.OO3r 2 c M , where 
r pc is the distance in parsec between the dust screen and the continuum source. 
The intrinsic dust model reproduces the flux rise observed around 660 A in key 
quasar spectra quite well. We present indirect evidence of a shallow continuum 
break near 670 A (18.5 eV), which would be intrinsic to the quasar continuum. 

Subject headings: galaxies: intergalactic medium — large-scale structure of universe 
— galaxies: active — radiative transfer — ultraviolet: general 
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1. Introduction 

The spectral energy distribution (SED) 
of active galactic nuclei (AGN) contains a 
significant feature in the optical-ultraviolet 
(UV) region, known as "the big blue 
bump". As for the emission lines in AGN 
spectra, it is generally believed that pho- 
toionization is the excitation mechanism 
of the emission lines superimposed to the 
continuum. Photoionization calculations 
that reproduce the AGN line ratios and 
equivalent widths favor an ionizing SED 
that peaks (in vF v ) in the extreme-UV 
(e.g. Mathews & Ferland 1987; Binette 
et al. 1988; Ferland et al. 1996; Korista, 
Ferland & Baldwin 1997). Satellite ob- 
servations of distant quasars, however, 
showed that the big blue bump peaks in 
vF v around > 1000 A (Kest) 1 (O'Brien et 
al. 1988; Zheng et al. 1997; Telfer et al. 
2002). This finding is best illustrated by 
the composite quasar SED constructed by 
Telfer et al. (2002) (hereafter TZ02), which 
was obtained by co-adding 332 HST-FOS 
archived spectra of 184 quasars with red- 
shifts between 0.33 and 3.6. The com- 
posite (reproduced in Fig. 15 of § 5) is 
characterized by a mean near-UV index 
a u of —0.69 {F v oc i/*"), steepening to 
~ -1.76 ± 0.12 in the far-UV. The fa- 
vored interpretation by TZ02, Zheng et al. 
(1997), and Shang et al. (2004) is that the 
observed continuum steepening is intrin- 
sic to quasars. Intriguingly, using archived 
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data from FUSE, whose sensitivity win- 
dow extends further in the UV, Scott et 
al. (2004) performed a similar compilation 
for 'nearby' active galactic nuclei (AGN) 
with redshifts z q < 0.7 and report the lack 
of any evidence of a steepening in the far- 
UV! Furthermore, the FUSE composite 
spectrum for nearby AGN is significantly 
harder than that of Telfer et al. (2002), 
with a„ ~ -0.561;! in the far-UV. Ar- 
guably, since nearby AGN are on average 
less luminous, they may possess an intrin- 
sically different SED. However, a detailed 
optical-UV study of a subset of the FUSE 
sample by Shang et al. (2004) fails to re- 
veal any correlation between the far-UV 
index and blackhole mass. Therefore, even 
though the HST-FOS and FUSE samples 
do not represent equivalent AGN popu- 
lations, the absence of steepening in the 
far-UV for the nearby sample cannot be ex- 
plained alone by difference in AGN popu- 
lations. A plausible explanation for the on- 
set and increasing importance of the break 
in distant AGN can be provided by inter- 
galactic absorption, since it would scale 
with distance. In an earlier paper, Binette 
et al. (2003) explored the possibility that 
the break might be the result of H I scat- 
tering by a tenuous intergalactic compo- 
nent that the authors associated with the 
warm-hot intergalactic medium. The mod- 
els, however, predicted a significant flux 
discontinuity in the region 1050-1190 A 
(A 6s.), which is not observed in quasar 
spectra, as shown by FUSE (e.g. Kriss et 
al. 2001). Furthermore, the warm-hot in- 
tergalactic medium is too ionized to pro- 
duce the amount of H I absorption needed 
to reproduce the break. In this Paper, 
we explore an alternative interpretation 
that is based on a different opacity vector, 
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namely dust, either intrinsic to the quasar 
environment or intergalactic. The vector 
responsible for the absorption will consist 
of grains made of carbon atoms, a ma- 
jor constituent of the interstellar medium 
dust, albeit here in crystalline form (nan- 
odiamonds). We will assume that the in- 
trinsic quasar SED consists of a simple 
powerlaw, and that deviations from the 
powerlaw are caused by absorption from 
crystalline carbon dust, either as pure cu- 
bic diamonds, or of the type observed in 
carbonaceous chondrite meteorites (e.g. 
Allende, Orgueil and Murchison) . Many 
mechanisms have been proposed to explain 
the formation of diamond nanocrystallites 
(c.f. §7). About half of them require 
intense UV irradiation. Interestingly, a 
significant UV flux is present in the two 
Herbig Ae/Be objects, for which nanodia- 
mond emission bands have been confirmed 
first (Van Kerckhoven et al. 2002). Pro- 
cesses that form nanodiamonds by UV ir- 
radiation are particularly relevant, since 
quasars are UV powerhouses and their en- 
vironment might lead to physical condi- 
tions that favor the emergence of carbon- 
based nanocrystallite grains. 

The paper is structured as follows: fol- 
lowing the introduction in § 1, we describe 
in § 2 the dust models and the algorithm 
used to compute the transmission function. 
The methodology and classification of the 
spectra are described in § 3. In § 4 and § 5, 
we present the intrinsic and intergalactic 
dust absorption models, respectively, and 
compare them with the observed spectra. 
In § 6 we decide, which of the two models, 
is to be preferred, and we discuss a possi- 
ble final model. In § 7 we focus on the for- 
mation and physical properties of nanodi- 



amonds, and follow with the conclusions in 
§8. 

2. Procedure and calculations 
2.1. Dust extinction curves 

In order to account for the sharp SED 
break by way of dust absorption, we looked 
for an absorption vector that peaks in the 
far-UV (A < 1000 A), and yet causes neg- 
ligible absorption at wavelengths longer 
than 1200 A. Ideally, as it is the case with 
the interstellar medium (ISM) dust, the 
grain particles should be composed of the 
most abundant elements. In both aspects, 
the crystalline form of carbon is the most 
appealing candidate and is the basis of this 
Paper. We will consider two types of mate- 
rials: the terrestrial cubic (pure) diamonds 
and the nanodiamonds as found in mete- 
orites. 

A comparison of the UV extinction 
properties of the terrestrial diamonds and 
the meteoritic nanodiamonds can be found 
in Mutschke et al. (2004). The authors 
separated the nanodiamonds from the Al- 
lende 2 meteorite sample and determined 
their optical constants. The meteoritic 
nanodiamonds differ in their optical prop- 
erties from the cubic diamond clS du result 
of chemical impurities (e.g. H, N) and of 
restructured or unsaturated bonds at their 
surface. 

Following a standard procedure, we cal- 
culated a set of dust extinction curves. We 

2 Carbonaceous chondrite meteorites are relatively 
rare, at a frequency of only ~ 3.5%. The Allende 
meteorite who fell on Earth near the town of Al- 
lende in the state of Chihuahua, Mexico, on Febru- 
ary 8th, 1969, is one of the most studied meteorites 
of its kind. 



3 



assumed that the grains are spherical and 
that a powerlaw describes the differential 
distribution of grain sizes, dn gr (a)/da = 
CgrUHcfii where n gr is the volume den- 
sity of grains, uh that of hydrogen, a the 
grain radius, ( the powerlaw 3 index and 
C gr the normalization constant, such that 
the density of grains becomes normalized 
to the abundance of the dust constituents 
with respect to hydrogen. We adopted 
the tabulated complex refraction indices 
n + ik of Mutschke et al. (2004) for the 
Allende meteorite nanodiamonds, and of 
Edwards & Philipp (1985) for the cubic 4 
(terrestrial) nanodiamonds. The Mie the- 
ory was used to compute the extinction 
cross-section Q ext (a, X,n,k), using a mod- 
ified version of the published subroutine 
BHMIE of Bohren & Huffman (1983). The 
extinction cross-section is normalized with 
respect to the gas density n#, using the 
following integrals: 

n H 0\ = n H C gr r ira 3+(; Q ex t(a, A, n, k) d\na 
V qr =|f na^dhia 

where a m j„ and a max are the minimum and 
maximum values of the grain radii con- 
sidered. The gas opacity is given by the 



integration of dr{ xt 



3 A log-normal distribution is more appropriate 
(Lewis et al. 1989) for describing the nanodiamond 
size distribution. However, because we intended 
to explore a much broader size range than that 
found in Allende, we considered it more convenient 
to use a powerlaw for this purpose. The reason is 
that only one parameter needs to be varied (a max ), 
while in the case of the log-normal distribution, we 
would need to vary two parameters simultaneously 
(the mean size value and the distribution's width). 

4 The extinction cross-section for cubic diamonds 
were extrapolated in the far-UV, since we could 
not find published laboratory values of the refrac- 
tion indices shortward of 413 A. 



/ A — n H a^dr. Neglect- 
ing the contribution of elements other than 
carbon to the composition of the nanodia- 
monds, we adopt a mean molecular weight 
of \x gr = 12 for the grains. The value of the 
normalization constant C gr is obtained by 
solving the following: 

Zcl^gr m H — PgrCg r Vg r 

where Zq< is the carbon abundance by num- 
ber with respect to H, m# the hydrogen 
atom mass and p gr the density of the grain 
material. The values adopted for p gr are 
2.3 (Lewis et al. 1989) and 3.51 gem -3 , for 
the Allende and the cubic nanodiamonds, 
respectively. For the grain size exponent, 
we adopted ( = —3.5. As we neither 
know the gas metallicity nor the dust-to- 
gas mass ratio, we assume all carbon is 
locked in dust and we adopt the solar value 
of 3.63 x 10~ 4 for the C abundance, for the 
sole purpose of procuring a convenient nor- 
malization. 

The continuous lines in Fig. 1, repre- 
sent the extinction curves adopted in this 
Paper. The curves labeled Dl and Al 
represent the "small size regime" extinc- 
tion curves for terrestrial diamonds (red 
curve) and for meteoritic nanodiamonds 
(blue curve), respectively. In both extinc- 
tion models, a m j n =3A and a max =25 A. In 
this regime, decreasing further a max would 
not alter the extinction curve. Our extinc- 
tion curve Al (blue line) is very similar 
in shape to the mass absorption coefficient 
curve determined by Mutschke et al. (2004) 
(see their Fig. 7 or our scaled version of it, 
the green dotted curve in Fig. 1). A third 
dust model, which is also useful, is the 
Allende curve labeled A3 (orange curve), 
whose grain sizes extend up to 200 A. The 
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peak cross-sections for the curves Dl, Al 
and A3, occur at wavelengths 640, 741 and 
787 A, respectively. Meteoritic nanodia- 
monds are known to possess a median ra- 
dius a of ~ 15 A. When increasing a max 
to a value of ~ 50 (75) A for the Allende 
(terrestrial cubic) nanodiamonds (respec- 
tively), one finds that the peak absorption 
starts shifting noticeably to the right and 
the absorption profile widens somewhat. 
This is illustrated by the two long-dashed 
line curves in Fig. 1, both calculated with 
a max =100 A. The above mentioned curve 
A3 further extends the grain size range to 
a max =200 A, which significantly shifts the 
broad absorption peak toward longer wave- 
lengths. 

As shown in §4, the extinction curves 
Dl and Al (or A3 in §5) can reproduce 
the wide range of continuum steepening 
observed in the far-UV in quasar SEDs. 
Both types can induce a sharp absorp- 
tion break, although the cubic diamond is 
more extreme in this respect. This does 
not occur with ISM dust extinction. For 
comparison, we plot an ISM dust model 
from Martin & Rouleau (1991) (with C = 
—3.5) in Fig. 1, which consists of silicate 
and graphite grains of sizes comprised be- 
tween a max = 2500 A and a min = 50 A 
(black short-dashed line). It is evident that 
the customary ISM extinction curve, while 
reaching a maximum in the UV, still ab- 
sorbs significantly longward of the peak, 
which gives rise to a shallow change of in- 
dex, rather than a sharp break. Grain size 
is not the main cause for such differences in 
relation to nanodiamonds, but rather the 
type of material being considered. To il- 
lustrate this, we show a small grains ex- 
tinction curve used by Binette, Magris & 



Martin (1993) to study the scattered con- 
tinuum of Pks2152— 69; it has the same 
composition as the ISM model, but the 
size range is reduced to a max = 500 A (black 
dotted line in Fig. 1). The cross-section 
redward of the extinction peak remains too 
shallow to reproduce the sharp break of 
quasars. Shang et al. (2004) explored the 
possibility that the standard ISM or even a 
SMC-like extinction curve could reproduce 
the QSO break. Their conclusions is that 
reddening by ISM or SMC-like grains "is 
not able to produce the spectral break seen 
in the AGN sample, without leaving a clear 
signature at longer wavelengths" (which 
usually is not seen). Notice that the cubic 
diamond curve shows a rather narrow peak 
at ~ 650 A followed by a lower plateau at 
~ 500 A. These particularities of the cu- 
bic diamond cross-section, together with 
the very steep rise shortward of 1000 A are 
unique features, which should leave a clear 
imprint, whenever this material is respon- 
sible for the extinction. The extinction at 
optical wavelengths due to nanodiamonds 
is negligible. For instance, with an opacity 
of unity at 912 A, the extinction in the V 
band (A rest ) is as small as Ay = 10~ 4 and 
5 x 10 _7 mag in the case of the Al and Dl 
curves, respectively. 

The extinction curves Dl and Al (or 
A3) as defined above will suffice to test the 
dust absorption hypothesis. Instead of us- 
ing optically known materials, one could 
have treated the absorption hypothesis as 
an inverse problem, working out the ex- 
tinction curve that succeeds best. Con- 
sidering that the current study is mostly 
exploratory in nature, we consider that it 
confers a higher degree of plausibility to 
use an empirical curve such as that of the 
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Allende meteorite, rather than an invented 
cross-section. 

2.2. Calculation of the transmission 
curve 

The basic assumption behind the cur- 
rent work is that the break observed in the 
spectra is a manifestation of dust absorp- 
tion and is therefore not an intrinsic fea- 
ture of the SED. A key aspect in evaluating 
how well the dust absorption hypothesis 
fares is to assume that we can extrapolate 
the powerlaw observed in the near-UV to 
the region underlying the break. Any de- 
parture of the observed spectrum from the 
extrapolated powerlaw will be modeled as 
"absorption". Only in §6.3 will a broken 
powerlaw be considered for the far-UV. In 
our notation, the "true" or intrinsic quasar 
SED is described by either one of the ex- 
pressions: 

F? = A(^) +a », or 

where z/ and Ao (=912 A) are the ioniza- 
tion thresholds of hydrogen in frequency 
and in wavelength units, respectively, 
while a u and f3\ are the corresponding 
powerlaw indices. A and B are normal- 
ization constants, one of which is set to 
unity according to whether F£ (i.e. A) or 
Fl (i.e. B) is plotted (respectively). In 
keeping with the tradition in AGN litera- 
ture, the index that we quote in the text 
will always be a u . In par with the work 
of TZ02 and Zheng et al. (1997), we pre- 
fer to plot F\(\) for most figures. The 
far-UV region beyond the break is better 
represented using F\ than F u , as can be 
appreciated by comparing panels a and b 
of Fig. 2. 



2.2.1. Intrinsic dust absorption 

The (modeled) transmitted flux in the 
one-dimension case is given by F™ od = 
T\F%, where T\ is the transmission func- 
tion, which for a point source is simply the 
exponential e~ r ^ . To compute the opac- 
ity r{ xt = N H a^ in the case of dust at the 
redshift of the quasar, all that is required is 
to specify the absorption column Nh and 
select one, or a combination, of the extinc- 
tion curves described in §2.1. 

2.2.2. Intergalactic dust and the simula- 
tion of the composite spectrum 

If the dust is intergalactic, it is neces- 
sary to integrate the transmission along 
the line-of-sight to the quasar. Because we 
also intend to simulate the process of con- 
structing a composite spectrum from syn- 
thetic quasar SEDs, we developed the fol- 
lowing numerical procedure. Briefly, the 
simulation of the composite will consist in 
multiplying each synthetic quasar spectra 
by the appropriate transmission function 
and then co-add them in the quasar rest- 
frame. The synthetic spectra, before dust 
absorption, share the exact same SED, but 
differ in redshift and in spectral coverage. 
In the simulation, we adopted the same set 
of quasar redshifts as those in the TZ02 
sample as well as the same set of wave- 
length limits, for the synthetic spectra, 
as those characterizing the TZ02 archived 
spectra. Each synthetic spectrum corre- 
sponding to a given quasar at redshift z q 
is divided into energy bins, and for each 
rest-frame bin Xj, we calculate the (mod- 
eled) transmitted flux F^ od = F«.T Xj = 
F^.e~ T ^ making use of the integrated 
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opacity along the line-of-sight up to z q : 

r(X 3 ) = f« ' n H (z) of(^) f dz 

where is the dust extinction cross- 
section evaluated at wavelength Xj/(l + z), 
and n H (z) the intergalactic dust density 
expressed in terms of the hydrogen den- 
sity. 

For the calculations of distances, ^ and 
baryonic densities, we assume the con- 
cordance ACDM cosmology with parame- 
ters derived from the WMAP experiment 
(Spergel et al. 2003), that is tt A = 0.73, 
tt M = 0.27, h = 0.71 with h = H /100 
and a baryonic mass of VL har h 2 = 0.0224 
corresponding to an hydrogen density at 
zero redshift of n° ar = 2.06 x 10~ 7 cm~ 3 . 

3. Methodology 

3.1. The initial database 

The spectral database adopted in this 
work is that of TZ02, which was kindly 
lent to us by R. C. Telfer. It comprises 332 
spectra, mostly HST-FOS 5 , of 184 quasars, 
already reduced and corrected for Galactic 
dust extinction. The spectra furthermore 
have been corrected by TZ02 for the pres- 
ence of Lyman limit absorbers (down to 
t > 0.3) as well as of the Lya absorption 
valley (caused by the cumulated absorp- 
tion from unresolved Lya forest lines). 

3.2. Near and far-UV spectral in- 
dices 

We define the far-UV as the wavelength 
region shortward of the break from 300- 
1000 A, the near-UV as the 1000 to 3200 A 

5 Note that the TZ02 sample includes 3 HST-STIS 
and 6 HST-GHRS spectra. 



region longward of the break, and the 
optical-UV as the 3200-4200 A region. 

Throughout this Paper, we will refer to 
the powerlaw index longward of the break 
as ajyuv, an d that shortward as apuv- We 
will assume that the intrinsic SED pow- 
erlaw index, a u , has the same value in 
the region of the break as in the near-UV, 
hence a u =a^uv- Whenever possible, the 
adopted value for a^uv wm be the value 
that we estimate empirically, using the ad- 
jacent near-UV region of the HST-FOS 
spectrum. This value is to be preferred 
over published values, which correspond to 
a forced fit of the combined optical-UV re- 
gion. The HST-FOS a NUV indices are usu- 
ally significantly harder. They are more 
appropriate for the exercise at hand, which 
relies on having a dependable SED descrip- 
tion immediately longward of the break 
that can be extrapolated one octave short- 
ward, in the region of the break itself. 

3.3. Pre-analysis of the HST-FOS 
sample 

Following a preliminary analysis of the 
TZ02 sample and of the properties of dust 
models, we established the following. (I)- 
Individual quasar spectra provide stronger 
constraints to the models than a single 
composite spectrum. The process of co- 
adding varied spectra to construct the 
composite inevitably lead to a loss of valu- 
able information. Modeling the compos- 
ite spectrum is probably an essential ex- 
ercise, but does not constitute a determi- 
nant proof of the validity of any model. 
For these reasons, we concentrate here on 
fitting individual quasar SED. (II)- Since 
combined multigrating spectra extend over 
a larger wavelength domain, they provide 
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stronger constraints for the models than 
single grating spectra. For this reason, 
this work considers only those 106 spec- 
tra of the TZ02 sample that correspond 
to a combination of two or three HST- 
FOS gratings 6 . (HI)- In the process of 
looking for patterns among the numerous 
spectral shapes encountered, we found it 
beneficial to classify these according to the 
signs by which dust absorption apparently 
manifests itself. The proposed classifica- 
tion is nothing more than a convenient and 
simplified characterization of the big blue 
bump phenomenology found among the 
archived HST-FOS spectra. By no means 
it implies that the quasars themselves are 
intrinsically different as a result of their 
spectra belonging to one class or another. 

3.4. Classification of multigrating 
spectra into classes A-D 

A physical insight on how dust can al- 
ter the continuum shape and account for 
the break has led to the classification of 
the multigrating spectra into four groups. 
The three most relevant groups are quali- 
tatively described in panel a of Fig. 2. The 
4 classes are defined as follows. 

(A) - The spectra that show a contin- 
uum steepening near 1000 A (A rest ) 
belong to class (A). The near-UV 
spectrum is hard in these spectra 
and the far-UV shows a moderately 

6 Of particular interest are quasars in the redshift 
range 0.9-2, for which the spectrum corresponds 
to a combination of 3 gratings. In these, the break 
is in full view and, in most cases, there is sufficient 
wavelength coverage, longward of the break, to in- 
fer the spectral index a^uv an d, shortward of the 
break, to distinguish absorption features that the 
models must reproduce. 



steepened continuum. PG 1148+549 
(z q = 0.969) can be considered the 
archetype of this class (see Fig. 5). 
We tentatively assign the 7 spec- 
tra (usually hight redshift quasars), 
whose near-UV FOS spectrum is not 
available longward of 1300A, to class 

(A) . HS 1700+6416 (z = 2.722) for 
instance is classified as class (A) 
(Fig. 13). More than 60% of quasars, 
whose spectra extended sufficiently 
into the far-UV to determine a^uv 
belong to class (A) alone. This may 
explain, why the spectra of this class 
individually resemble the TZ02 com- 
posite shape, since the composite is 
after all the result of averaging spec- 
tra that most often than not belong 
to class (A). 

(B) - The spectra that show a sharp break 

near 1000 A, followed by an ex- 
tremely steep continuum drop short- 
ward of the break, belong to class 

(B) . The near-UV spectrum is hard 
in these spectra. PG 1248+401 (z q = 
1.03) in Fig. 3 can be considered the 
archetype of this class. Another ex- 
ample is Pks 0122-00 (z q = 1.07) in 
Fig. 4. Objects in this class are not 
that common (only 6), though strik- 
ing 7 by their lack of a significant flux 
in the far-UV. 

(C) - The spectra that show a continuum 

that is already soft longward of the 
break, that is, up to > 1600 A, be- 
long to class (C). The soft region 
of the continuum now extends to 



7 Ton 34, which was reported to have an index of 
&fuv = —5.3 by TZ02, is another example but 
only a single grating spectrum exists. 
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include the continuum beneath the 
CivAA1549 doublet (or even down 
to Cin]AA1909 in some cases). A 
representative class (C) spectrum 
is the quasar 1130+106Y in Fig. 11 
(z q = 0.54). In many cases, a single 
powerlaw does not fit the near-UV 
part well and, in other cases, the 
index is very steep (< —1) through- 
out the whole spectrum, as exem- 
plified by 3C279 in Fig. 11 (green 
spectrum). In the far-UV, these ob- 
jects show characteristics of either 
class (A) or (B), that is, they are ei- 
ther flat in F x or very steeply declin- 
ing, as illustrated by MC 1146+111 
(z q = 0.863, blue spectrum), which 
exhibits a class (B)-like break. We 
found 8 objects with the above char- 
acteristics 8 . 

(D) - The high redshift quasars (3 objects) 
that we could not make sense of, 
belong to class (D). They are ob- 
jects that show an inflection or wide 
through in the far-UV. HE 1122-1649 
is one example (see blue spectrum in 
Fig. 18). We do not rule out that the 
troughs could be associated in some 
cases to one or more Lja absorption 
systems. 

Only 61 of the available 106 multigrating 
spectra extended sufficiently shortward of 
the break, that is, down to at least 900 A, 
to ensure proper classification. Therefore, 

'Having access only to the near-UV spectrum may 
suffice for classifying quasars into class (C) but, in 
keeping with the above classification rules, we did 
not consider nor classify any spectrum that did 
not extend down to at least 900 A in the far-UV. 



only this subset of 61 quasars has been an- 
alyzed in detail and modeled. Of these, 44 
are class (A), 6 class (B), 8 class (C) and 
3 class (D). We will mainly focus on class 
(A) and (B) spectra. These two groups to- 
gether represent 82% of the classified ob- 
jects and will suffice for the purpose of test- 
ing the dust absorption hypothesis. Many 
interesting spectra that could not be shown 
in this Paper will appear elsewhere (e.g. Bi- 
nette et al. 2005a, b, c). As for class (C), 
dust appears to be related to some of the 
observed characteristics of at least a frac- 
tion of them (§ 4.3) but further work will be 
needed to reach definite conclusions. The 
few objects that form class (D) are puz- 
zling and will not be modeled with dust 
absorption in this Paper. 

3.5. Absorption models considered: 
intergalactic vs intrinsic 

In order to explore how dust absorption 
might be the real cause of the observed 
break, a decision must be made on where 
the dust is located. The answer to this 
question defines two basic 9 types of ab- 
sorption models: (a) the dust is intrinsic 
to the environment of the quasars, and (b) 
the dust fills the intergalactic space. In 
case (a), the transmission function is de- 
rived directly from the extinction curve in 
the rest-frame of the quasar, as mentioned 
in §2.2.1. With this category of models, 
we may reasonably expect the amount of 
dust to vary more or less at random from 
object to object. In case (b), the dust dis- 
tribution is intergalactic and fills large vol- 
umes of space. We therefore expect that 

9 As the spectra have already been corrected for 
Galactic reddening, there is no reason to be con- 
cerned by Galactic dust absorption. 
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such a dust distribution should be, to a 
first order, homogeneous, since the dust 
becomes a cosmo logical component unre- 
lated to the quasars. With case (b), the 
models predict the same transmission for 
objects of comparable redshifts, indepen- 
dently of class. The dust density, uh(z), 
function of redshift as mentioned in § 2.2.2, 
has to be determined, requiring extra con- 
straints. Intergalactic models imply enor- 
mous amounts of dust, since it is cosmolog- 
ical. We will first study the intrinsic dust 
case and then proceed to the intergalactic 
case. 

In all figures, when overlaying a dust ab- 
sorbed model to the rest-frame spectrum 
of a particular quasar at redshift z q , we 
follow the following coding: the continu- 
ous line part depicts the wavelength region 
corresponding to an idealized FOS spec- 
trograph window extending from 1250 to 
3600 A (X bs.), while the dashed line rep- 
resents an extension into the far-UV down 
to 915A (X bs.), as would be available using 
the FUSE satellite. A dotted line is used 
outside these two observer-frame windows. 

4. The case for intrinsic dust 

The case in favor of intrinsic dust ab- 
sorption is best made by going through 
each class in order of increasing complexity 
of the dust model that it requires, that is 
in order B, A and C. 

4.1. Class B spectra 

Although class (B) objects are not nu- 
merous, they gave us an important clue on 
how to disentangle various effects result- 
ing from dust absorption. What charac- 
terize this class is the very steep drop of 



the UV flux shortward of 1000 A (A rest ). 
Class (B) objects can easily be accounted 
for by simply using the extinction curve 
Dl consisting of terrestrial cubic nanodi- 
amonds and adjusting as needed the ab- 
sorption column 10 N H . This is illustrated 
in Fig. 3, which shows the spectrum of the 
archetype class (B) quasar, PG 1248+401. 
The red line corresponds to a model us- 
ing the curve Dl and an absorption col- 
umn density N H of 3.2 x 10 20 cm~ 2 (here- 
after, the notation A*2o= 3.2 will be used). 
The assumed underlying power law index 
is a u = a^uv — 0.0, which is the index 
that best fits the emission-line-free contin- 
uum longward of Lja. In all our plots, it 
is the quasar spectrum that we scale, until 
an overlap with the model is obtained in 
the near-UV. The resulting spectrum scal- 
ing factor, Mf 14, is listed in each caption 
in units of 10 14 erg -1 cm 2 s A. Since we as- 
sume the intrinsic SED to be described by 
a simple powerlaw (until § 6.3), the plotted 
models in all figures will correspond to the 
function F A mod = T A x (A/912)-( 2+c H and 
the y-axis can be used to infer the trans- 
mission value T\ for any value of A. 

The only free parameter of the above 
Dl dust model is the column A^o- The 
position of the peak in transmission over- 
laps surprisingly well with that of the ob- 

'It is likely that only a small fraction of carbon 
is actually locked into nanodiamond grains. Sup- 
posing we independently knew the dust-to-gas ra- 
tio due to nanodiamonds A^ TG and that it was 
smaller than 0.0031, which corresponds to deplet- 
ing all the carbon onto dust, then the absorption 
columns quoted in this paper would have to be 
multiplied by the factor 0.0031/A£ TG . The galac- 
tic ISM dust is characterized by a much larger 
Ajf^ ~ 0.009, since it contains many other 
atomic species than C. 
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served spectrum. This is the result of the 
very sharp drop in cross-section of curve 
Dl longward of ~ 1000 A (see red curve 
in Fig. 1), a property unique to cubic di- 
amonds among dust grains composed of 
pure carbon. It is important to note that 
class (B) objects cannot be accounted for 
by making use of the Allende extinction 
curves Al or A3, because these are charac- 
terized by a broader absorption peak and 
the peak itself is shifted toward higher A 
values. The blue line in Fig. 3 illustrates 
the case of using curve Al. The relative 
success of the Dl dust model is telling us 
that if dust were indeed responsible for 
the break in class (B) objects, it must 
mostly consist of cubic diamonds. This 
does not rule out that a small fraction of 
the dust may be of the Allende type. This 
is demonstrated by the green line, which 
corresponds to a model with A r 2 o=2.8 and 
a linear dust mixture of 85% of Dl grains 
and 15% of Al grains. Hereafter, we will 
use the notation of f m = 0.85 to repre- 
sent the fraction of Dl grains, with 1 — f D1 
being the fraction of Al grains. 

Other class (B) quasars are Pks 1229—02 
and Pks 1424—11, which are quite simi- 
lar to PG 1248+401. They require models 
with dust columns of N 2 o=3.6 and 2.0, 
respectively, and an extinction curve con- 
sisting totally or mostly of cubic diamond 
type Dl. In Fig. 4, we show another (B)- 
type spectrum, quasar Pks 0122— 00, which 
requires somewhat less dust. For this ob- 
ject, the green line model has iV 2 o=2.0 and 
consists of a dust mixture with /di=0.8. 
It appears that the model with a mixture 
of dust types Dl and Al results in a supe- 
rior fit than the blue line model with pure 
cubic diamond extinction (with / D1 =1.0 



and iV 2 o=2.3). 

In summary, class (B) spectra require 
dust that predominantly consists of cubic 
diamonds. The particular absorption char- 
acteristics of cubic diamonds fit the ob- 
served far-UV steep drop particularly well. 
Due to the large dust opacities implied, 
there subsists little or no far-UV flux to 
be observed in these objects shortward of 
800 A. If photoionization by high energy 
UV photons is the excitation mechanism of 
the emission lines, it is puzzling to find that 
the same high excitation emission lines are 
observed in class (B) quasars, devoid of 
hard UV, as in other objects that do have a 
hard continuum (e.g. HS 1700+6416 with 
OiFuv=— 0.55; Reimers et al. 1989). A pos- 
sibility is that the dust lies outside the 
Broad Line Region (BLR). In this case, 
the emission line BLR clouds would be ex- 
posed to an ionizing continuum that is not 
absorbed. 

4.2. Class A spectra 

4-2.1. A mixture of the two nanodiamond 
grain types 

Class (A) spectra show a continuum 
break that is far less pronounced. The 
far-UV continuum is often flat 11 in F\, 
with an index cxfuv of order —1.7 short- 
ward of 1000 A. One may reasonably ex- 
pect that the absorption dust columns are 
simply smaller than in the previous case. 
This is confirmed by models. Unlike class 
(B) objects, where one species of dust is 
clearly favored, class (A) objects gener- 
ally require an extinction curve that corn- 
By 'flat' we mean an approximately horizontal 
continuum segment in an F\ plot, which trans- 
lates into an index j3\ ~ 0, that is, a. v ~ —2. 
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bines the extinction from Allende nanodi- 
amonds with that of cubic diamonds, that 
is, nanodiamonds with and without sur- 
face impurities (§2.1). The spectra of 
PG 1148+549 in Fig. 5 will serve to illus- 
trate this point. All models shown have 
the same column of A^ 20 =1.05 and the 
same SED with a v —— 0.2. Either extinc- 
tion curve Dl or Al can give rise to a 
flat continuum immediately longward of 
the break, but the onset of the break turns 
out inappropriate for both the Dl extinc- 
tion curve (red line) and the Allende Al 
curve (blue line), as seen in Fig. 5. An ex- 
tinction curve consisting of a mixture with 
f D1 =0.6 (60% Dl and 40% of Al grains) 
on the other hand provides quite an ac- 
ceptable fit to the break (green curve). We 
found that, by using a proper mixture of 
the two grain types, one can fit all class (A) 
objects. Three prominent emission lines 
stand out above the far-UV continuum in 
PG 1148+549, shortward of 900 A (the line 
identifications shown in the figures follow 
those proposed by TZ02). It is interesting 
to note that extinction by pure cubic dia- 
monds (red line) gives rise to a narrow dip 
near 700 A, a feature not observed in this 
particular quasar. 

4-2.2. Applying the intrinsic dust model 
to class (A) spectra 

A mixture of the two nanodiamond 
types is very successful in reproducing the 
break in all class (A) spectra. More specif- 
ically, we can fit all observed 1000 A breaks 
assuming one of the four following values 
of fm'- 0, 0.3, 0.6, and 1.0. A finer sub- 
division in most cases is not warranted by 
the data, because the spectra very rarely 
extend sufficiently in the far-UV for the fit 



to be sensitive to small changes infm- In 
a few cases, only the onset of the break 
is seen, and we could not determine with 
certainty whether the value of 0.3 or 0.6 is 
more appropriate. At any rate, the value 
°f /di— 0.3 appears to be the most frequent 
as indicated in the histogram of Fig. 6, but 
there remains a substantial fraction of ob- 
jects that require a different dust mixture. 

In Fig. 7, we present the distribution of 
gas columns derived from fitting class (A) 
and class (B) spectra. The mean N 2 o value 
for class (A) is 1.02 with a standard devi- 
ation of 0.29. There are two high-redshift 
spectra, for which there is no evidence of 
dust, with an upper limit of 0.1, which are 
HS 1700+6416 and HE 2347-4342 (they 
have not been included in the average). 
If we combine class (A) and (B) spectra 
into a single group and assume that they 
are part of the same population, we de- 
rive A^2o=l-20 and a standard deviation 
of 0.60. The distribution shows that the 
presence of nanodiamond dust is the rule 
in quasars rather than the exception. In 
many quasars, the amount of dust inferred 
is comparable. For instance, 39 quasars 
have a column 0.6 < N 2 o < 14, which rep- 
resents 78% of the 50 class (A+B) spectra. 

As for the distribution of a u describing 
the near-UV continuum, we obtain a mean 
value for class (A) of cxnuv = —0.44 with 
a dispersion of 0.21. This average consid- 
ers only the 22 objects, for which a reli- 
able estimate of cinuv could be determined 
directly from the HST spectra. It is sig- 
nificantly harder than the mean value of 
—0.69 reported by TZ02 and the median 
value of -0.83 ± 0.04 for local AGN re- 
ported by Scott et al. (2004), presumably, 
because the softer class (C) spectra are not 
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included in our average. 

4-2.3. Dust models predict a rise short- 
ward of ~ 700 A 

Due to the rapid decrease in the dust ex- 
tinction cross section in the far-UV, short- 
ward of the cross section peak (see Fig. 1), 
an inescapable feature of dust absorption is 
that a rise in transmitted flux always oc- 
curs in the far-UV, shortward of ~ 700 A. 
Of the 7 multigrating spectra that ex- 
tended down to 600 A and showed evidence 
of dust absorption, we found evidence of a 
sharp flux rise in 4 of them. This test was 
not conclusive for the remaining 3 spectra. 
The observation of a steep rise in the far 
UV is the strongest evidence in favor of 
dust absorption and will be presented in 
more detail. 

The first spectrum with a sharp rise is 
PG 1008+1319, which is shown in Fig. 8. 
We adopt the value of +0.13 as near-UV 
index, as determined by Neugebauer et al. 
(1987) . The green line corresponds to an 
intrinsic dust model with /di=0.3 while 
the gray line corresponds to /di=0.6. The 
column in both models is A^2o=l-2. Clearly 
the model with an extinction dominated 
by Allende nanodiamonds (red line) gives 
a better fit to the break. Reducing f^i 
further would cause the break 's onset to 
occur at too long a wavelength (e.g., the 
/di=0 model in Fig. 5). 

A second example, is Pks0232— 04 of 
Fig. 9. A fit of the near-UV continuum fa- 
vors a v indices in the range —0.2 to —0.4. 
To be definite, we adopt the steeper SED 
with —0.4. We verified that the same con- 
clusions are reached when using the harder 
index. The red line model, which is more 
satisfactory, corresponds to pure Dl dust 



with/ D1 =1.0 and A^o=0.90 while the green 
line model corresponds to /di=0.8 and 
-^20=0.93. Due to the predominance of cu- 
bic diamonds, a narrow dip at 650 A stands 
out in models with f m > 0.7. This dip ap- 
pears to be saddled by two prominent emis- 
sion lines, O ill and NeVlll, both of which 
are also visible in the composite SED of 
TZ02, but not as prominently. Interest- 
ingly, Scott et al. (2004) discuss the na- 
ture of a narrow dip seen blueward of the 
NeVlll emission in their near-AGN com- 
posite. The interpretation they favor is 
that of blueshifted absorption by NeVlll. 
Another explanation might be that absorp- 
tion by cubic nanodiamonds is responsible 
for this feature. Even though the Dl dip is 
broader, it might partially be filled by the 
NeVlll emission line. 

A third example is provided by the 
much higher redshift quasar, HS 1307+4617, 
at z — 2.129, which is plotted in Fig. 10. 
There is no HST-FOS spectrum that covers 
the near-UV. Instead, we adopt the value 
of aNuv=0.0 as inferred from a spectrum 
of D. Reimers and reproduced in Koratkar 
& Blaes (1999). The three models super- 
imposed to the spectrum in Fig. 10 have 
the same column A^ 20 =1.3 and differ only 
by their proportion of the Dl and Al dust, 
as follows : gray line foi—0.8, green line 
/di=0.6, and purple line /di=0.3. The 
green line model with /di=0.6 provides a 
better fit. It is interesting to note that the 
disjoint part (yellow segment) of the GHRS 
spectrum (taken with grating G140L) is 
not consistent with the far-UV extrapo- 
lation of the models. A solution to this 
problem is found in §6.3. A fourth ex- 
ample is 1623.7+268B, shown as the dark 
green spectrum at the bottom of Fig. 9. 
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To summarize, the intrinsic dust model 
is phenomenologically very successful, as 
it can not only reproduce the break within 
the dominant class (A) spectra, but it can 
also account for the continuum sudden rise 
in the far-UV, in the 650-700 A region. 
The extinction curve that is required to 
model the far-UV consists of a mixture of 
the two nanodiamonds grain types: the cu- 
bic diamonds and the Allende type. 

4.3. Class C spectra 

This class is defined by the near-UV 
shape, which is very soft, taking the ap- 
pearance of a flat continuum in F\ or, 
in some cases, of a bump longward of 
1200 A. We have already presented (§3.4) 
the three examples shown in Fig. 11. The 
frequent rounded shape of the spectrum 
from near to far-UV suggests the possi- 
bility that, in some cases at least, the 
spectra might be reddened by dust simi- 
lar to that of the Galactic ISM. One such 
case is 1 130+1 06 Y (black line spectrum in 
Fig. 11), which we tentatively model using 
the ISM dust model of Martin & Rouleau 
(1991) (see black dashed line extinction 
curve in Fig. 1). It consists of a mixture of 
silicate and graphite grains. The magenta 
line is a better model in the far-UV. It 
combines ISM extinction (80%) with that 
of cubic diamonds (20% of Dl grains by 
mass). The column is A^" 2 o=9.0 and the 
assumed SED has a v =— 0.25. Note the 
significant attenuation of the continuum, 
which reaches a factor of 5 near 1000 A. 
The optical extinction in the V band due 
to ISM-like dust is relatively modest, with 
Ay = 0.4 mag. Not including the con- 
tribution from nanodiamonds extinction 
would result in a flatter continuum without 



much of a break, as illustrated by the yel- 
low line model (pure ISM extinction with 
N 20 =9.0). 

Evidence of reddening by ISM dust ap- 
pears to be present in class (C) spectra, but 
nanodiamond dust is nevertheless required 
to explain the break when it is present. A 
few class (C) spectra show a steep drop 
shortward of 1000 A as in class (B) spectra. 
We have not explored the possibility of a 
combination of the three extinction curves 
Al, Dl and ISM. It is possible that some 
class (C) spectra possess a SED (from near 
to far-UV) that is intrinsically steeper than 
in other classes, with a u in the range —2 
to -1.2. 

5. The case for intergalactic absorp- 
tion 

The fact that the absorption columns 
take on similar values in the intrinsic case 
for the numerous class (A) spectra, invites 
us to explore the hypothesis that the dust 
pervades the intergalactic space instead of 
being confined to the environment of each 
quasar. Following this hypothesis, the dis- 
tribution of the dust bears no relation to 
the quasars, but is a function of distance 
(i.e. z). By the same token, we expect 
the grains composition to be more uniform 
in the intergalactic case than in the intrin- 
sic case. The intergalactic model does not, 
however, imply that there cannot be an ad- 
ditional intrinsic dust component local to 
some quasars, as appears to be the case for 
class (B) and (C) spectra. On the other 
hand, the case for intergalactic dust will 
be more convincing if a minority of class 
(A) spectra require additional absorption 
above the one provided by the intergalac- 
tic model. 
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The predictive value of the intergalactic 
model resides in the function chosen to de- 
scribe the dust density with redshift. Such 
a function will not only allow the model- 
ing of the break in individual quasars but 
also the simulation of the composite SED 
derived by TZ02, following the procedure 
defined in §2.2.2. 

5.1. Constraining the dust behavior 
with redshift 

As a working hypothesis, let us assume 
that the dust is intergalactic and for now 
consider only class (A) spectra. Since the 
absorption occurs along the line-of-sight to 
each quasar, its impact can extend over 
the whole far-UV domain as a result of 
the redshift effect and cosmological expan- 
sion. In a Universe that evolves and ex- 
pands, any cosmological quantity such as 
the density of the Lja absorbers, the den- 
sity of of quasars, the star forming rate, 
etc, is known to evolve strongly with red- 
shift, that is with time. The same must ap- 
ply to the hypothesized intergalactic dust. 
There should exist an epoch (z p ) at which 
the dust density reaches a peak. To de- 
scribe such a peak, we adopt a paramet- 
ric form for the dust density 12 n^(^), sim- 
ilar to that used by Baldry & Glazebrook 
(2003) to describe the cosmological star 
formation rate. It consists of a broken pow- 
erlaw joining at redshift z p : 

, < _ n%(l + z) e for z < z } 

Uh[ ~ z >- n° H (l + z p y-i{l + Z y forz>z } 

where n° H is the density at zero redshift, e 
and 7 are the low and high redshift indices 

The density nn as defined in this work is not a 
co-moving but rather a local quantity. 



and z p the intersection of the two truncated 
powerlaws. Hereafter, we will use n° in 
units of 10~ 8 cm~ 3 to express the density 
at zero redshift. 

To constrain the parameters describing 
the function tih(z), we proceeded as fol- 
lows. Since Scott et al. (2004) did not find 
evidence of a continuum break in nearby 
AGN (z q < 0.7), this suggests that the 
peak in absorption occurred at an earlier 
epoch rather than in the local Universe. 
A positive index for e, in which absorp- 
tion increases with look-back time, will 
have the effect of reducing the importance 
of the 1000 A break within the local Uni- 
verse. Another indication of the increase 
in the importance of the break with red- 
shift can be appreciated in Fig. 12, where 
we plot OLpjjy as measured for each quasar 
by TZ02. To derive the mean values repre- 
sented by the squares, we distributed the 
measured indices into 5 redshift bins and 
then calculated the average aNuv within 
each bin. (The continuous line simply con- 
nects the 5 mean values.) After trial and 
error and varying e, we found that similar 
fits to the break could be obtained, using 
any value within the interval 1.5 < e < 3.5. 
To constrain e more effectively would have 
required including the nearby AGN obser- 
vations with the FUSE satellite (Scott et 
al. 2004). To be definite, we hereafter 
adopt the value e=+2. 

To determine the behavior of iih(z) 
at the other redshift end, we can com- 
pare the very high redshift quasar spectra 
(z > 2.5) with those at intermediate red- 
shifts. As it turns out, most high- z HST- 
FOS spectra are single-grating and can't 
be used for that purpose. Fortunately, 
there exist two quasars with high S/N 
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and wide spectral coverage that we could 
analyze in greater detail, HS 1700+6416 
and HE 2347-4342, which are plotted as 
F v in Fig. 13 13 . The powerlaw indices for 
HS 1700+6416 and HE 2347-4342 that we 
adopt are a u = —0.55 (from Reimers et al. 
1989) and +1.70 (inferred from the best 
model), respectively. The missing param- 
eter values defining uh{z) were arrived at 
using various constraints, as described be- 
low. 

Because of the redshift effect, the HST- 
FOS spectra of both HS 1700+6416 and 
HE 2347-4342 do not cover the typical 
break region at 1000 A (A rest ). Viewed 
from the perspective of intergalactic dust, 
however, the absorption break should have 
shifted to shorter wavelengths (with z), as 
demonstrated below in §6.2. As a con- 
sequence, the continuum shape's depar- 
ture from that of a pure powerlaw must 
be the result of the hypothesized inter- 
galactic dust, if such a model is to be of 
any use. Despite the ragged appearance 
of both continua in Fig. 13, caused by the 
many absorption systems along the line- 
of-sight, it is clear that both show a gen- 
eral curvature 14 or change of index, which 
intergalactic absorption must be able to 
explain. It turns out that such a cur- 
vature can be reproduced by intergalactic 
dust models using either of the extinction 
curves, Al or Dl. To select the appropri- 
ate extinction, we required that the dust 

Because of the hardness of HE 2347-4342 ((3\ ~ 
—3), using F° bs is much more convenient than 

pobs 

Note that if we attempt to fit the above curvature 
using intrinsic dust, the absorption actually goes 
the wrong way, making the transmitted spectrum 
appear even harder, as illustrated by the green line 
model calculated with 7V 2 o=0.5 and /di=0.3. 



model successfully reproduced the break 
observed in the lower redshift spectrum 
of PG 1148+549, which is the archetype 
of class (A). This second constraint effec- 
tively rules out cubic diamonds as shown 
by the red line model of PG 1148+549 in 
Fig. 14. In the case of intergalactic mod- 
els with dust curve Al, the break occurs 
at somewhat too short a wavelength (see 
blue line in Fig. 14). To compensate for 
the redshift smearing effect, we extended 
the grain size range of the Allende nanodi- 
amonds, extending it up to a max =200 A. 
This defines the new extinction curve A3 
(see orange curve in Fig. 1) used in all our 
intergalactic calculations. This extinction 
curve A3 was found to provide an overall 
better fit to class (A) objects, in the inter- 
galactic case. 

Having selected the optimal extinction 
curve, strong constraints on z p and n° H 
can now be derived by varying these pa- 
rameters until an acceptable fit of the two 
high-z quasars is found. We found that 
7 is loosely constrained to negative values 
< —1.4. To further constrain the function 
tih(z) and 7, we made use of the compos- 
ite quasar SED of TZ02 (shown in Fig. 15). 
We simulated this composite by co-adding 
synthetic dust-absorbed SEDs of the same 
redshifts and spectral widths as those in 
the actual TZ02 sample, following the pro- 
cedure described in §2.2.2. This exercise 
indicated a preference for somewhat larger 
values of 7 than the range favored above 
using the two high-z quasars. To be def- 
inite, we adopted the value of 7 = —1.5, 
which corresponds to the overlap between 
the two types of constraints. The result 
of combining these different constraints in 
an iterative fashion has been that an ac- 
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ceptable fit to the broad curvature of both 
high- ,2 quasars of Fig. 13 occurs when us- 
ing the value z p ~ 0.4 for the peak dust 
redshift. The resulting orange line model, 
which requires ng=3.4, now overlays both 
continua in Fig. 13 as well as the outline 
of the break in PG 1148+549 (Fig. 14). In 
conclusion, the intergalactic dust model 
can account for the progressive steepening 
of the powerlaw index observed shortward 
of 500 A in HE 2347-4342. 

5.2. Applying the intergalactic model 

A reasonable expectation of the inter- 
galactic dust model is that it should apply 
to all the classes defined in §3.4. This is 
not to say that additional absorption by 
a local dust component cannot take place 
in some quasars. For instance, class (B) 
quasars, although dominated by intrinsic 
dust, as shown in §4.1, can also be mod- 
eled as the sum of intergalactic absorption 
and absorption by intrinsic Dl dust. An 
example of such complementarity is given 
by class (B) Pks 0122-00 in Fig. 4. An or- 
ange line model is plotted, which includes 
absorption by both intergalactic A3 dust 
and intrinsic Dl dust of column A^ 2 o=l-0. 
Except toward the far-UV, this orange line 
model (mostly hidden by the foreground 
red line!) is almost identical to the pre- 
viously described green line model (§4.1), 
which consisted of an intrinsic dust mix- 
ture with /di=0.8 and column A r 20 =2.0. 
A similar comparison can be established 
with the intergalactic orange line model of 
PG 1248+401 (Fig. 3). 

When attempting to simulate the com- 
posite SED of Fig. 15, it turns out that n§ 
must be increased, from 3.4 to 4.7. If not, 
the model lied significantly above the com- 



posite. Such a model with n§ increased to 
4.7 is shown by the purple continuous line 
in Fig. 15. This simulation though imper- 
fect is encouraging, if we consider that our 
simulation assumed a single powerlaw in- 
dex a v =— 0.6, while the TZ02 composite 
sampled widely different energy distribu- 
tions. Furthermore, TZ02 have combined 
spectra of classes (A)-(D), while the pure 
intergalactic model is intended for class 
(A) objects. It is presumably for that rea- 
son that a significant improvement is ob- 
tained, as shown by the orange line model 
in Fig. 15, when one combines intergalac- 
tic absorption with intrinsic absorption by 
a column as small as A" 20 =0.15 of Dl dust. 

We have repeated the same exercise as 
for the intrinsic case (§4.2) of applying the 
intergalactic model to all class (A) spectra. 
We found that 20 spectra could be reason- 
ably fitted with the proposed value n§=3.4 
(the standard model), while 14 spectra re- 
quired either an increase of ~ 10% in n® 
or the addition of an intrinsic column of 
dust, most frequently of type Dl. Further- 
more, 4 and 6 spectra required n§ to be 
reduced (approximately) to 2.4 and 3.0, re- 
spectively. The database is therefore not 
entirely consistent with the expectation of 
an homogeneous distribution of intergalac- 
tic dust. 

6. Merits of intergalactic vs intrin- 
sic models 

Which model is to be preferred? We will 
compare the merits and problems of each 
type of model and, in conclusion of this 
section, we will present a final model for 
the 1000 A break of quasars. 
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6.1. The intrinsic dust hypothesis 

We have shown that intrinsic dust mod- 
els could account not only for the break, 
but also for the flux rise at shorter wave- 
lengths (e.g. Figs. 8-10). Overall, the in- 
trinsic model is extremely successful across 
the whole class (A) and class (B) samples. 
The fit to the far-UV rise in HS 1307+4617, 
on the other hand, is not entirely satisfac- 
tory. 

Is it possible to simulate the TZ02 com- 
posite assuming only intrinsic dust? One 
difficulty is that the intrinsic model is mute 
about how other parameters like N 2Q , fm 
or a v might vary with increasing z. How- 
ever, since each class (A) spectrum can be 
fitted quite well by varying f D1 and N 20 , 
which is an approach that we consider su- 
perior to that of simply fitting the TZ02 
composite, it can then be argued that not 
being successful in the simulation of the 
composite is of secondary importance. Al- 
though we consider this to be true, we nev- 
ertheless attempted to simulate the com- 
posite, because it reveals real trends in 
quasar SEDs. The fact that the TZ02 com- 
posite remains very soft at very short wave- 
lengths instead of showing a far-UV rise, 
must be explained somehow. 

The silver dashed line in Fig. 15 illus- 
trates our initial attempt to simulate 15 
the composite SED, assuming a u =— 0.6 
and keeping all the input parameters con- 
stant with z. The column is N 2 o=l-0 and 
/di=0.3. The far-UV flux is obviously pre- 

'In a situation in which none of the parameters 
defining the intrinsic model varies with redshift 
(as for the silver dashed line model in Fig. 15), all 
the absorbed SEDs are identical and there is no 
need to co-add the spectra in order to simulate 
the composite. 



dieted too strong. This happens, because 
the extinction cross-section falls off too 
rapidly at very short wavelengths, and a 
steep rise in F\ becomes unavoidable short- 
ward of 600 A. The simulated composite 
simply tends toward the slope given by the 
index a u . If we vary the column with red- 
shift, by defining a function N 2 q(z), we ob- 
tain the absurd result that, in order for the 
simulated composite to overlap the TZ02 
composite, the dust column would have 
to increase sharply with redshift. This 
is not only an ad hoc dust behavior, it is 
also contradictory to the absence of any in- 
crease with z q of the columns determined 
in §4.2. In addition, it is at odds with the 
lack of absorption in the two high-z spec- 
tra of HS 1700+6416 and HE 2347-4342, 
for which we determined absorption upper 
limits of A^2o < 0.1. Attempts to model 
the curvature in these two spectra with in- 
trinsic dust result in absorption features at 
the wrong end of the spectrum. In effect, 
local dust makes both spectra appear even 
harder than they already are, as illustrated 
by the green line model in Fig. 13 calcu- 
lated with iV"2o=0.5 and / D1 =0.3. 

6.2. The intergalactic dust hypoth- 
esis 

The few shortcomings mentioned above 
for the intrinsic case disappear with inter- 
galactic dust. By construction, the con- 
tinuum of the two high-redshift quasars, 
HS 1700+6416 and HE 2347-4342, can be 
reproduced. On the other hand, the break 
can be fitted only for a qualified majority 
of class (A) spectra, while the other spec- 
tra usually requiring intrinsic dust to be 
added to the model. The TZ02 composite 
can be reproduced, albeit with a density 
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n° increased by 40% (for which we have no 
satisfactory explanation to propose). The 
far-UV index olfw-, when evaluated at the 
fixed wavelength of 800 A, exhibits the cor- 
rect trend with redshift, as shown by the 
long-dashed line in Fig. 12. One may argue 
that the amount of dust implied by the in- 
tergalactic model is excessive if not plainly 
unreasonable, but it is not an impossible 
amount. The fraction of the baryonic mass 
that the value rig=3.4 corresponds to is 
17% (see §6.3) , assuming that the mean 
cosmic carbon metallicity is about solar 
and that the dust is intergalactic, because 
it was expelled from galaxies by radiation 
pressure (Ferrara et al. 1991) or through 
super novae of type II. 

The intergalactic model, on the other 
hand, makes stringent predictions about 
how the break ought to shift (and soften) 
with increasing redshift. This is shown 
in Fig. 16, in which the transmission func- 
tion is plotted at representative z q values. 
The continuous part of each T\ curve cor- 
responds to the fiducial spectrograph win- 
dow of 1250-3600 A (A o6s .) [see §3.5] and 
shows what part of the break is visible at 
a given redshift z q . Notice that when the 
redshift exceeds values of ~ 1.5, the break 
is markedly shifted toward shorter wave- 
lengths. Of the three spectra presented in 
§4.2.3, which showed a clear flux rise in 
the far-UV, only one is of sufficiently high 
redshift to test this, HS 1307+4617 with 
z = 2.129. It's spectrum is shown again in 
Fig. 17 and can be compared with the pure 
intergalactic model with n°=3.4, which is 
represented by the brown dashed line. The 
SED is the same as earlier, that is a u =0.0. 

The gradual break (or curvature) seen 
in the brown dashed line model in Fig. 17 



not only occurs at very short wavelengths 
(~ 550 A), but is extremely shallow. Obvi- 
ously, in order to fit the sharp break char- 
acterizing the HS 1307+4617 spectrum, 
additional intrinsic absorption must be 
considered. Such a model is represented 
by the orange line, which is a model that 
combines intergalactic with intrinsic dust. 
The local dust column is A^o=0.8 with 
a dust composition /di=0.6. The fit to 
the observed flux rise is surprisingly good, 
much better even than with the pure in- 
trinsic case represented by the green line 
in the previous Fig. 10 of the same quasar. 
Even more suggestive is the disjoint spec- 
trum obtained with GHRS grating G140L 
(yellow spectrum in Fig. 17), which de- 
spite its lower S/N appears to prolong the 
far-UV rise of the multigrating spectrum 
(black line). The intergalactic model is 
marginally consistent with the continuum 
level set by this spectrum segment, in con- 
trast with the pure intrinsic model, which 
rises too steeply (see green line in Fig. 10). 
The green dashed line in Fig. 17 repre- 
sents the contribution of intrinsic dust ab- 
sorption that is present in the orange line 
mixed model. In summary, even though in- 
trinsic dust is the main contributor to the 
sharp break observed in HS 1307+4617, 
the signature at the shortest wavelengths 
expected in high-z spectra result of 
intergalactic dust appears to be indepen- 
dently confirmed in this quasar. 

Is the intergalactic hypothesis vindi- 
cated? It turns out not to be the case. 
In effect, a rather poor fit is provided 
for the two other quasar spectra that 
showed a far-UV flux rise, PG 1008+1319 
and Pks 0232-04. As is the case for 
HS 1307+4617, these two quasars require 
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additional intrinsic absorption, with dust 
columns (A^o) of 0.5 and 0.25 and dust 
mixtures (/di) of 0.0 and 1.0, respectively. 
However, even when combining this addi- 
tional absorption with intergalactic dust, 
the far-UV rise cannot be reproduced at 
all, as shown by the corresponding orange 
lines in the Fig. 8 of PG 1008+1319 and in 
the new Fig. 18 of Pks 0232-04. The far- 
UV continuum level is predicted too low in 
the Pks 0232-04 model and, in both fig- 
ures, the flux rise (orange lines) occurs at 
too short a wavelength. 

6.3. Evidence of a higher energy 
break? 

As indicated above, the discrepancy 
of the intergalactic model for reproduc- 
ing the far-UV rise in PG 1008+1319 and 
Pks 0232— 04 could not be resolved. This 
inadequacy of the model is sufficiently sig- 
nificant to reject the intergalactic dust hy- 
pothesis at the assumed density 16 . More- 
over, in more than one aspect, the inter- 
galactic model is implausible. Assuming 
that the mean metallicity of galactic mat- 
ter (stars and interstellar gas) at current 
epochs is near solar, as derived by Calura 
& Matteucci (2004a, b), then the fraction 
of cosmic carbon required is n° H /n® ar = 
3.4 x 10~ 8 /2.06 x 10~ 7 = 0.17. Apart from 
the unreasonable fraction of cosmic carbon 
that must exist in crystalline form (17%), 
it would require an improbable fine tuning 

ie If 

we adopt an intergalactic model that uses signif- 
icantly less dust than ng=3.4, it would be at the 
cost of having more intrinsic dust present and this 
in a larger fraction of class (A) spectra, if not the 
majority of them. Hence, such a model would not 
contribute in an essential way in explaining the 
break and would have to be discarded on account 
of Ockham's razor principle. 



so that this dust would not be accompa- 
nied by larger amounts of the more com- 
mon flavors (silicates, graphites, PAHs, ...) 
as one would expect if it was formed in 
supernovae and later expelled in the inter- 
galactic space. The analysis of cosmologi- 
cal supernovae has ruled out the existence 
of large quantities 17 of intergalactic dust 
of 'normal' composition (Perlmutter et al. 
1999). 

Ruling out the intergalactic model in fa- 
vor of the intrinsic model has an additional 
interesting consequence. The curvature 
present in the spectra of HS 1700+6416 
and HE 2347— 4342 must now be consid- 
ered an intrinsic feature of the energy dis- 
tribution in these two quasars, rather than 
the manifestation of intergalactic absorp- 
tion. 

We are left with the intrinsic dust hy- 
pothesis to account for the 1000 A break 
in quasars. What could be missing from 
this model so that the few remaining prob- 
lems it has could be resolved? We will 
hypothesize that the shallow rollover ob- 
served in the far-UV in HS 1700+6416 and 
HE 2347—4342 is a manifestation of a uni- 
versal high energy cut-off present in all 
quasar SEDs. To test this, we first model 
the cut-off in a way that does not depend 
on the a^uv of the underlying SED. This 
is achieved by singling out the transmission 
curve calculated at z q = 2.8 (correspond- 

We have calculated the extinction that ISM dust 
would produce if it followed the same intergalac- 
tic distribution as the nanodiamonds (§5.1) with 
rig=3.4, that is in amounts that correspond to 17% 
of the Galactic dust-to-gas ratio. We find that for 
an object at z = 0.5, the selective extinction is 
Eb-v = 0.022 (X bs.), much in excess of the value 
of 0.002 inferred from cosmological supernovae by 
Perlmutter et al. (1999). 



20 



ing to the average redshift value of the two 
high-z quasars) and considering it to be 
a valid description of the continuum soft- 
ening that might also apply to the other 
quasars. We obtained a parametric fit of 
this transmission curve at that redshift, us- 
ing the following expression: 

c x = (1 + [\/\ brk } fs y f 1 (2) 

where 5 is the power law index change and 
/ a form factor, that may vary from ob- 
ject to object. In essence, when a quasar 
SED is multiplied by Cx and 5 < 0, a shal- 
low steepening takes place at \b r k, which 
increments the underlying powerlaw index 
by 5. The sharpness of the cut-off is set 
by the form factor / > 1. Hereafter, we 
will assume that the true intrinsic quasar 
continuum is given by C\ x F%, where F% 
is the quasar powerlaw that we have been 
using so far and a v the index set by the 
near-UV continuum. To illustrate the role 
of the form factor, let us take the spec- 
trum (blue spectrum in Fig. 18) of class 
(D) quasar HE 1122— 1649 as an example, 
given that it shows a clear far-UV cut-off. 
Longward of the cut-off, the powerlaw in- 
dex is a v = —0.6, reaching —2.6 shortward 
of the break situated at ~ 670 A. Hence the 
parametric function C\ is characterized by 
5 = —2.0 and A& r fc=670A. The gray and 
magenta lines in Fig. 18 correspond to form 
factors of / = 2.8 and 10, respectively. 

Turning our attention to the T\ curve 
at a redshift z g =2.8 in Fig. 16, the para- 
metric fit using Eqn. 2 yields the values 
5 = -1.6, / = 2.8 and A ferfc =670A. It is 
shown by the thick gray dashed line. This 
cut-off can now be applied, not only to 
HS 1700+6416 and HE 2347-4342, but to 
any other quasar as well. It will result in a 



new SED that is characterized by the same 
near-UV powerlaw as before, but with a 
shallow far-UV cut-off centered on 670 A, 
that is, at 18.5 eV. 

If we incorporate the above far-UV 
rollover into our previous quasar SED of 
index a v = —0.6, we obtain the distribu- 
tion represented by the black dotted line 
in Fig. 19. Assuming this modified SED, a 
remarkable improvement in the simulated 
composite is now obtained, as shown by 
the dark gray continuous line in Fig. 19. 
(The silver dashed line model represents 
the earlier case without the far-UV intrin- 
sic cut-off). Further improvements of the 
simulated composite can be obtained by 
varying some of the input parameters with 
redshift. The study of these more complex 
models, however, exceed the scope of this 
paper and will be reported elsewhere. 

It is important to emphasize that the 
intrinsic model with a shallow cut-off fits 
markedly better two of the three quasar 
spectra that showed a far-UV rise. For in- 
stance, the earlier problem encountered 
(§4.2.3) in fitting the far-UV rise in 
HS 1307+4617 (green line in Fig. 10) has 
now disappeared, as shown by the cyan 
line model in Fig. 17, which assumes the 
SED with the above far-UV rollover. An 
improvement of the fit to the far-UV rise 
in Pks0232— 04 is also obtained, as shown 
by the cyan line model in Fig. 18. In the 
case of PG 1008+1319, it makes little dif- 
ference whether the high energy cut-off is 
there or not (compare the cyan and green 
line models in Fig. 8). 

In summary, the whole database as rep- 
resented by the TZ02 composite, as well 
as the individual spectra showing the far- 
UV rise, are both consistent with the pres- 
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ence of an 18.5 eV cut-off in quasars. This 
proposed new cut-off is being masked in 
the vast majority of quasars with redshift 
< 2.5 by the more prominent 1000 A break, 
which we believe is entirely due to nanodi- 
amond dust absorption. 

7. Nanodiamonds: the infrared-UV 
connection 

Nanodiamonds are, to date, the most 
abundant presolar grains, both in mass 
and numbers, that have been extracted 
from primitive carbonaceous meteorites 
(Mutschke et al. 2004, and references 
therein), but their detection in the ISM has 
been elusive. Diamond crystallite emis- 
sion bands in the 3.3-3.6 /zm region due 
to surface C-H stretching modes of hy- 
drogenated nanodiamonds have been es- 
tablished with confidence for a few Herbig 
Ae/Be objects and one carbon-rich post- 
AGB star HR4049 (Guillois et al. 1999; 
Van Kerckhoven et al. 2002; Acke & van 
den Ancker 2004). Van Kerckhoven et al. 
(2002) presented a detailed analysis of the 
ISO-SWS spectra of the two Herbig ob- 
jects, HD 97048, and EliasI, as well as of 
the post-AGB HR4049. They applied a 
physical model to the emission profile of 
the 3.53 fim band and inferred a temper- 
ature of 950 K and 1000 K for HD 97048 
and EliasI, respectively. Assuming radia- 
tive equilibrium between photoheating and 
far-infrared cooling for the grains, the au- 
thors could estimate the UV radiation flux 
impinging the diamonds in these three ob- 
jects. The diameter range they inferred for 
the crystallite diamonds is 2a ~ 10-100 A. 
Interestingly, the multiwavelength data for 
both HD 97048, and EliasI, as well as for 
the post-AGB star HR4049, indicate that 



the 3.53 /iin emission takes place within a 
disk-like structure. The distance between 
the star and the emission region is < 9 
and < 22 AU, in HD 97048 and EliasI, re- 
spectively. The formation site that these 
authors favor for the crystallite diamonds 
is in situ formation within the disk rather 
than within the ISM or via ejection from 
stars. In the field of AGN, sub-arcsecond 
VLT observations by Rouan et al. (2004) 
using NAOS+CONICA revealed wavelike 
structures in the mid-infrared, which the 
authors propose might be due to emission 
by nanodiamonds at a temperature close 
to sublimation. 

To explain the predominance of nan- 
odiamond grains in primitive meteorites, 
several formation mechanisms have been 
proposed, such as: (a) chemical vapor de- 
position from stellar outflows (Lewis et 
al. 1987), (b) impact shock metamorphism 
driven by supernovae (Tielens et al. 1987), 
(c) energetic ion bombardment by a super- 
nova (Daulton et al. 1995), (d) UV an- 
nealing of carbonaceous grains (Nuth & 
Allen 1992), (e) nucleation in organic ice 
mixtures by UV photolysis (Kouchi et al. 
2005), and (f) chemical conversion of PAH 
clusters to nanodiamonds in the presence 
of UV radiation (Duley & Grishko 2001). 
It is interesting to note that the last three 
processes involve UV radiation. The above 
post-AGB and the two Herbig Ae/Be stars 
emit UV radiation (Van Kerckhoven et al. 
2002), a fact which is possibly related to 
the formation of the observed nanodia- 
monds. 

Could a similar formation process op- 
erate within the UV-intense environment 
of quasars? The indication that cubic 
diamonds dominate in class (B) quasars 
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might be related to an evolutionary se- 
quence of the grains. A possible scenario 
might be the following. Via the process 
of dehydrogenation of PAH clusters by 
quasar UV radiation, hydrogenated nan- 
odiamonds form, with optical properties 
similar to the Allende type. If the UV radi- 
ation heats up the nanodiamonds beyond 
1300 K, a process of surface dehydrogena- 
tion begins, which may cause the grain op- 
tical properties to become more similar to 
that of cubic diamonds. Finally, the disap- 
pearance of H-stretch cooling may result in 
a runaway heating, followed by graphitiza- 
tion and eventually to sublimation of the 
grains. 

To confirm the existence of nanodia- 
mond grains in AGN, one could attempt 
to detect the far-infrared emission bands 
caused by hydrogenated nanodiamonds 
(Van Kerckhoven et al. 2002; Jones et al. 
2004a). However, AGN are intrinsically 
very strong far-infrared emitters, and the 
signature of any narrow emission band will 
certainly be diluted. For instance, the dust 
silicate feature at 9.7/zm predicted by cal- 
culations (Laor & Draine 1993) is not ob- 
served as often as expected in AGN. The 
UV radiation absorbed by nanodiamonds 
represent at most 10% of the energy inte- 
grated over the whole SED. Let us assume 
the AGN unification picture with a bi-cone 
opening angle of 45°. If the nanoparticles 
are located outside the BLR, then a frac- 
tion of only 0.10 x (1 - cos ) < 0.01 
of the quasar bolometric luminosity will 
be reprocessed into far-infrared emission 
by nanodiamonds. Assuming a uniform 
covering factor of unity within the radia- 
tion bi-cone, the minimum dust mass re- 
quired by the intrinsic model is given by 



0.044A" 2 o Tp C x (1 - cos ^-) M , where r pc is 
the distance in parsecs separating the dust 
screen from the central UV source. For an 
arbitrary distance of one parsec and N 2 o of 
unity, the dust mass implied is 0.0033 M Q , 
a value independent of the assumed dust- 
to-gas ratio (see footnote 10). 

To the extent that Allende-type nan- 
odiamonds might be a candidate carrier 
(see Jones et al. 2004b) for the Extended 
Red Emission observed in nebulae between 
5400 and 9500 A, it is conceivable that a 
fraction of the UV flux absorbed by dust 
might be re-emitted by photoluminescence. 

8. Conclusions 

We have presented evidence that in- 
dicates that dust absorption by nanodi- 
amonds is successful in reproducing the 
1000 A break as well as the far-UV rise 
seen at shorter wavelengths. Could the 
agreement between the intrinsic dust mod- 
els and the spectra be simply the result 
of a coincidence between the break loca- 
tion and the far-UV extinction proper- 
ties of nanodiamonds? To rule out such 
a possibility will require that an inde- 
pendent confirmation of the presence of 
nanodiamonds be found. Detection of 
grain emission in the far- infrared, at 3.43 
and 3.53 fim, is one possibility, although 
this emission mechanism works only with 
surface-hydrogenated grains (see §7). An- 
other possible route would be to observe 
selected quasars in order to extend the UV 
coverage in objects for which only the on- 
set of the break is seen so far. The idea 
would be to look for a confirmation of a 
flux rise shortward of 700 A in as many 
quasars as possible. This would require 
high quality observations using a satellite 
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with far-UV sensitivity. One possibility 
might be the R=1000 spectrometer on- 
board the projected World Space Observa- 
tory satellite, which is expected to offer a 
sensitivity window covering the range If 00 
to 3500A (X obs ) (Barstow et al. 2003). 

We have ruled out that the dust caus- 
ing the f 000 A break is predominantly in- 
tergalactic on the basis that it is not re- 
quired per se and that the far-UV rise could 
not be modeled satisfactorily using inter- 
galactic dust. Furthermore, the amount of 
crystalline carbon that is needed turns out 
to be improbably large. Since intergalac- 
tic dust is not responsible for the contin- 
uum rollover observed in HS 1700+6416, 
HE 2347-4342, at 2 ~ 2.8, we have pro- 
posed that this feature is a manifestation 
of a higher energy break, near 18.5 eV, 
which is presumably intrinsic. Including 
the same break in the other quasars SED 
markedly improve the simulation of the 
composite as well as the detailed modeling 
of the far-UV rise in PG 1008+1319 and 
Pks0232— 04. At any rate, such a break 
is bound to take place somewhere in the 
far-UV so that the quasar SED connects 
smoothly with the soft-X rays. In effect, 
the optical- X-ray index, aox, which re- 
lates the monochromatic continuum flux at 
2500 A to that at 2 keV, is characterized by 
values in the range 1.3 to 1.6 (equivalent to 
an a v between —1.3 and —1.6). Given that 
our mean a nuv index for class (A) is much 
harder, with —0.44, a continuum turnover 
must take place somewhere in the far-UV. 
Our results suggest that such turnover oc- 
curs at 18.5 eV. We are currently in the 
process of studying how the far-UV and 
the soft X-rays may join together (Haro- 
Corzo et al. 2005). 



Intrinsic dust models require gas columns 
of the order 10 20 cm -2 , assuming solar C 
abundance and full depletion onto nanodi- 
amond grains. For comparison, in the so- 
lar neighborhood, a U-band extinction of 
a tenth of a magnitude by ISM dust corre- 
sponds to a gas column of 1.9 x 10 20 cm~ 2 
(Whittet 2002). Larger columns (but same 
dust masses) would be implied for our 
models, if we assumed smaller dust-to-gas 
ratios. To the extent that the nanodia- 
mond dust lies in the vicinity of the AGN, 
the minimum dust mass required by in- 
trinsic models is small ~ 0.003^0^ M 
(see §7). This value is independent of the 
assumed dust-to-gas ratio, because the ab- 
sorption gas columns scale inversely with 
it (see footnote 10). 

We have shown evidence that ISM-like 
dust might be playing a role in explaining 
the continuum appearance of a fraction of 
quasars, the so-called class (C) quasars. To 
confirm this suggestion will require more 
complex models than presented here, in 
which more extinction components might 
have to be contemplated (ISM, SMC, nan- 
odiamonds, SiC, ...). An alternative is that 
class (C) quasars emit with an intrinsically 
much softer SED. 

How does fare the intrinsic dust hypoth- 
esis in relation to the result of Scott et al. 
(2004), in which the break is apparently 
absent in local Universe AGN? An interest- 
ing scenario comes to mind. An inspection 
of the 45 multigrating spectra, which could 
not be classified, because their spectra did 
not extend down to 900 A (A rest , see §3.4), 
reveals that the softer class (C) spectra are 
more frequent at lower redshifts. Interest- 
ingly, a fraction of class (C) quasars do not 
show any break (see spectrum of 3C279 in 



24 



Fig. 11). A possibility might be that there 
is a secular evolution of the dust proper- 
ties, with nanodiamonds being absent in 
quasars with z q > 2.5, later becoming pre- 
dominant at z q < 2, and, finally, being 
progressively replaced by ISM-like dust for 
z q < 0-7. 

Interpreting the 1000 A break in terms 
of dust absorption may contribute to re- 
solve the following issues in the AGN field: 

(a) - The continuum rise in the 650- 

700 A region, seen in a few individual 
quasar spectra, is not predicted by 
any accretion disk model or any other 
continuum emission model known to 
the authors. Such a rise on the other 
hand is expected with nanodiamond 
dust absorption. 

(b) - If the 1000 A break is intrinsic to the 

ionizing continuum of quasars, the 
mean ionizing photon energy turns 
out rather small, making it difficult 
to account for the observed luminosi- 
ties of the high excitation lines (e.g. 
Korista, Ferland & Baldwin 1997). 
With the alternative interpretation 
of dust absorption, the SED turnover 
is pushed to higher energies and the 
break is an artifact of line-of-sight 
dust absorption. 

(c) - The puzzling fact that the high exci- 

tation emission lines in UV deficient 
quasars (e.g. class (B) quasars such 
as PG 1248+401 in Fig. 3) are com- 
paratively as luminous as in other 
quasars. This is easily explained by 
dust absorption provided the dust 
lies outside the BLR or under any ge- 
ometry, in which dust only affects the 



observer's line-of-sight and not the 
BLR line-of-sights to the UV source. 

(d) - The far-UV rise observed in a few 

quasars, HS 1307+4617, PG 1008+1319 
and Pks0232— 04, must be followed 
by a SED turnover at higher ener- 
gies. A detailed modeling of this 
flux rise, assuming dust absorption, 
as well as the shallow rollover seen 
in HS 1700+6416 and HE 2347-4342 
(and the sharper cut-off in class 
(D) HE 1122-1649, Fig. 18) are both 
consistent with the presence of a con- 
tinuum cut-off at 18.5 eV. Even if 
such a cut-off is not directly percep- 
tible in individual spectra in the rest 
of the sample, it is consistent with 
the far-UV slope seen in the TZ02 
composite. 

(e) - The narrow continuum dip short- 

ward of the NeVlll emission line in 
the composite spectra of TZ02 and 
Scott et al. (2004) is not accom- 
panied by a similar absorption dip 
blueward of O VI in emission, as one 
might expect if the trough was due to 
absorption by an outflowing ionized 
wind. Instead, the NeVlll trough 
could be the result of the narrow 
absorption peak that characterizes 
the cubic diamonds (Dl) extinction 
curve. 
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Fig. 1. — Extinction cross-sections for nan- 
odiamonds of radii in the range 3-25 A 
from the Allende meteorite (blue contin- 
uous line labeled Al) and for dust grains 
consisting of (terrestrial) cubic diamonds 
(red continuous line labeled Dl). The two 
long-dashed curves illustrate the effect of 
increasing a max to 100 A, corresponding to 
dust models A2 and D2 (labels not shown). 
Finally, the curve A3 corresponds to the 
case of increasing a max to 200 A (contin- 
uous orange line). The dotted section of 
the terrestrial diamond curves shortward 
of 413 A corresponds to an extrapolation, 
as the refraction indices are not available. 
The green dotted curve is the Mutschke 
et al. (2004) mass absorption coefficient 
curve, which is renormalized so that its 
maximum coincides with the Al curve. 
The two curves barely differ, except long- 
ward of 1200 A. The black dashed curve 
corresponds to a model of the ISM dust 
by Martin & Rouleau (1991). The "small 
grains" dotted curve is the same model, 
but with a max reduced (from 2500 A) to 
500 A (Binette, Magris & Martin 1993). 

Fig. 2. — A qualitative description of the 
three main spectral classes defined in § 3.4. 
Panel a: arbitrary scale in F u , panel b: 
arbitrary scale in F\ for the same spec- 
tra shown in panel a. The F\ represen- 
tation offers a clearer view of the steep- 
ened far-UV region. The dotted line il- 
lustrates typical variations within a given 
class. Of the 106 multigrating HST spec- 
tra available, only the 61 spectra that ex- 
tended beyond the break, down to at least 
900 A, could be classified reliably. The ver- 
tical dashed line represents the position of 
the Lyman limit (912 A X res t). 



Fig. 3.— Spectrum of PG 1248+401 
(Xrest)- The combined spectrum (grat- 
ings G190H+G270H) of PG 1248+401 is 
shown by the continuous thin black line 
and has been multiplied by the scaling fac- 
tor 0.75 x 10 14 erg" 1 cm 2 s A (hereafter, the 
notation M/i 4 =0.75 will be used). Pointers 
indicate the position of relevant emission 
lines identified by TZ02 in their compos- 
ite spectrum. Red line: absorption model 
as a function of rest-frame A using extinc- 
tion curve Dl (cubic diamond), column 
N 20 =3.2 and a SED with a„=0.0. The in- 
trinsic quasar SED assumed (not shown) 
is a powerlaw F^ = (A/912)~( 2+Q ^), hence 
Fg 12 = 1 (in all figures). The blue line 
represents an absorption model of equal 
column ^20=3.2, which uses the extinc- 
tion curve Al (Allende nanodiamonds). 
The green line corresponds to a model 
with A" 2 o=2.8 and a dust mixture with 
/ m =0.85, that is 85% Dl grains and 15% 
Al grains. Underneath the green line lies 
the intergalactic model (orange line) of 
§5.1 combined with intrinsic dust with col- 
umn ^20=1-8 and f D1 =1.0. In all mod- 
els, the continuous line part corresponds 
to a fiducial spectrograph window extend- 
ing from 1250 to 3600 A (X obs .), while the 
dashed line represents an extension of the 
model into the far-UV, down to 920 A 
(X bs.)- A dotted line is used in models out- 
side these two observer-frame windows (see 
§3.5). 
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Fig. 4.- Spectrum of Pks 0122-00 
(gratings G190H+G270H) multiplied by 
Mf i4— 1.1. The notation is the same as in 
Fig. 3. Red line: absorption model using 
extinction curve Dl (cubic diamond), col- 
umn A^2o=2.3 and an SED with a u =— 0.55. 
The blue line represents a model using 
the extinction curve Al (Allende nanodi- 
amonds) and a column iV"2o=2.0, while the 
green line corresponds to a model with 
iV2o=2.0 and a dust mixture /di=0.8. Un- 
derneath the green line lies the intergalac- 
tic model (orange line) of §5.1 combined 
with intrinsic dust with column 7V 2 o=l-0 
and / m =1.0. 



Fig. 7. — Histogram of the distribution 
among the 44 class (A) and 6 class (B) 
spectra of the gas columns, N20, results in 
an optimal fit of the UV break. The first 
bin at N 2 o=0-l corresponds to objects for 
which only an upper limit of N H could be 
determined. 



Fig. 5.— Spectrum of PG 1148+549 (grat- 
ings G130H+G190H+G270H) multiplied 
by Mj 14 =0.46. The notation is the same 
as in Fig. 3. Red line: intrinsic absorp- 
tion model using extinction curve Dl (cu- 
bic diamond) with a column A^o=l-05 
and a SED with a u =-0.2. The blue line 
represents the same model but assuming 
the extinction curve Al (Allende nanodia- 
monds), while the green line corresponds 
to a dust mixture model with f m =0.6. 
This mixed dust model provides a sat- 
isfactory fit of the continuum underlying 
the three emission lines: OlII, NeVIII and 
Oiii. 

Fig. 6. — Histogram of the distribution 
among the 44 class (A) spectra of the frac- 
tional contribution, fci, of cubic diamond 
dust (Dl) to the extinction curve, which 
results in an optimal fit of the UV break. 
The contribution of the meteoritic Allende 
nanodiamonds (Al) is 1 — fm- 



Fig. 8.— Spectrum of quasar 

PG 1008+1319 (gratings G150L+G270H) 
multiplied by Mj 14 =0.65. The notation is 
the same as in Fig. 3. Notice the far-UV 
flux recovery shortward of 720 A. The 
assumed near-UV index is a^=0.13 for all 
models (Neugebauer et al. 1987). Green 
line: absorption model with dust mixture 
/di=0.3 and a column iV"2o=l-2. The gray 
line represents a different dust mixture of 
/di=0.6 with same column iV2o=l-2 (it 
lies in the far-UV behind the green line) 
The orange line is the intergalactic model 
introduced in §5.1 to which intrinsic dust 
with /di=0.0, and column V 2 o=0.50 has 
been added. The cyan line is a model that 
assumes an SED, modified by a shallow 
cut-off C\ as defined in §6.3 (Eqn. 2), a 
dust mixture with /di=0.3 and a column 

#20=1.2. 
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Fig. 9.— Spectrum of Pks 0232-04 
(gratings G150L+G270H) multiplied by 
Mf 14 =0.73 (black line). The notation is 
the same as in Fig. 3. Notice the far-UV 
flux recovery shortward of 640 A. The as- 
sumed near-UV index is a„=— 0.4 for all 
models. Red line: absorption model of 
Pks 0232-04 with column iV 20 =0.90 and 
extinction curve Dl (/di = 1.0). The green 
line represents a different dust mixture 
with f m =0.8 and column iV 20 =0.93. The 
dark green line spectrum corresponds to 
quasar 1623.7+268B (M / i 4 =l.l), which 
also shows a rise in the far-UV. 



Fig. 11. — Three class (C) spectra: 
cyan line MC 1146+111 (M fu =1.3), green 
line 3C279 (M /14 =0.70) and black line 
1130+106Y (M /14 =0.46). The notation is 
the same as in Fig. 3. Both MC 1146+111 
and 3C279 have been suitably scaled so as 
not to overlap with the quasar 1130+106Y, 
which is being modeled. 3C279 is flat in 
F x , hence a u ~ —2. Yellow line: dust 
absorption model of 1130+106Y assum- 
ing a v =— 0.25 and an extinction curve cor- 
responding to pure Galactic ISM extinc- 
tion with iy 2 o=9.0. (This last model has 
been multiplied by 1.3 before plotting.) 
The magenta line represents a dust mix- 
ture of ISM-type grains (80%) and terres- 
trial diamond grains Dl (20%). The V- 
band extinction implied by this model is 
Ay = 0.4 mag. 



Fig. 10.— Spectrum of HS 1307+4617 
(gratings G190H+G270H, black line) mul- 
tiplied by M /14 =0.75. The yellow dis- 
joint part corresponds to a GHRS spec- 
trum with grating G140L (same M/14). 
The notation is the same as in Fig. 3. No- 
tice the far-UV flux recovery shortward of 
700 A. The near-UV index is 0.0 in all mod- 
els. Three models are shown that used the 
same column A^2o=l-3, but different dust 
mixtures: green line: /zn=0.6, the purple 
line: /di=0.3, and the gray line: /£>i=0.8. 



Fig. 12. — Far-UV powerlaw indices ctFuv 
of the quasar sample, as determined by 
TZ02, as a function of redshift. The solid 
squares (connected by a continuous line) 
represent an average of a^jjy within 5 red- 
shift interval bins. The long-dashed line is 
the spectral index at the fixed wavelength 
of 800 A, calculated for the standard inter- 
galactic model (n§=3.4). 
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Fig. 13.— Spectrum in F v of 

HS 1700+6416 (panel a, M /14 =0.75) and 
HE 2347-4342 (panel b, M /14 =7.8) as a 
function of A. The yellow disjoint segment 
corresponds to a GHRS spectrum with 
grating G140L (M /14 =7.8). Orange line: 
absorption model assuming extinction 
curve A3 and intergalactic distribution of 
dust with n§=3.4 (in units of 10~ 8 cm" 3 ), 
z p =0A, e=+2 and 7=— 1.5 (see §5.1). 
The quasar energy distribution assumed 
is described by the magenta dotted lines 
and correspond to a u =—0.55 (Reimers et 
al. 1989) and +1.70, for HS 1700+6416 
and HE 2347-4342, respectively. The 
green dashed lines illustrates the effect of 
having intrinsic rather than intergalactic 
dust. The column is iV2o=0.5 and the 
dust composition is /di=0.3. Increasing 
this column results in selectively more 
absorption at the longer wavelength end, 
the opposite of what is required. 



Fig. 14.— Spectrum of PG 1148+549 mul- 
tiplied by Mj 14 =0.46. The notation is the 
same as in Fig. 3. The models represent 
intergalactic absorption calculations, using 
three different extinction curves and as- 
suming the same SED with a u —— 0.2 as in 
the earlier Fig. 5 of the same quasar. Blue 
line: extinction curve Al, red line: extinc- 
tion curve Dl, and orange line: extinc- 
tion curve A3. The best fit is provided by 
the orange line model based on the curve 
A3, for which the grain sizes extend up to 
a max =200 A (see Fig. 1). The values of z v 
that provide an acceptable fit to the two 
high-redshift quasars of Fig. 13 are 0.8, 0.6 
and 0.4 in the cases of the models that use 
the curves Dl, Al and A3, respectively. In 
these models, e=+2 and 7= — 1.5, as de- 
fined in § 5.1. 
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Fig. 15. — Black thin line: composite en- 
ergy distribution derived by Telfer et al. 
(2002) and scaled here by a factor 0.61. It 
combines both radio-loud and radio-quiet 
quasars. Purple line: simulated composite 
absorption model assuming the dust is ex- 
clusively interg atactic with n§=4.7 (§5.2). 
Each synthetic spectrum, before being ab- 
sorbed and averaged, had the same SED 
consisting of a powerlaw of index a u =— 0.6, 
which is represented by the magenta dot- 
ted line. The spectral window limits 
(A & s .) and the corresponding redshift of 
each spectrum, before being averaged, was 
taken from the sample of TZ02. The or- 
ange line is the same intergalactic model 
but combined with absorption by a col- 
umn, A^2o=0.15, of dust intrinsic to each 
quasar (/zn=1.0). Silver dashed line: sim- 
ulated composite absorption model assum- 
ing the dust is intrinsic and not inter- 
galactic. The column is A^ 20 =1.0 and the 
mixture /di=0.3, both remaining invariant 
with redshift. 



Fig. 16. — Transmission functions T\(z q ) 
for z q values of 0.1, 0.2, 0.5, 1.0, 1.5, 2.0, 
2.8,3.5 and 5.0, assuming the standard in- 
tergalactic dust model and the extinction 
curve A3. The parameters used to de- 
fine the function nn{z) are z p =0A, e=+2 
and 7= — 1.5 (Eqn. 1 in §5.1). The con- 
tinuous part of each curve corresponds to 
the fiducial spectrograph window of 1250- 
3600 A (X bs.) at the corresponding quasar 
redshift. The thick gray dashed line is a 
fit of the z q = 2.8 transmission curve using 
the cut-off function C\ defined in Eqn. 2, 
with A brfe =670A, <5=— 1.6 and / = 2.8. 



Fig. 17. — A repetition of the spectrum 
of HS 1307+4617 of Fig. 10. In all mod- 
els, the assumed near-UV index is 0.0, as 
in Fig. 10. The brown dashed line corre- 
sponds to the standard intergalactic model 
introduced in §5.1. The orange line is 
the same intergalactic model, but com- 
bined with an intrinsic dust column of 
A^2o=0.8 with /di=0.6. The green dashed 
line shows the contribution by intrinsic 
dust in the case of this combined inter- 
galactic model (orange line). The cyan line 
is the intrinsic dust model, but assuming 
that the powerlaw is modified by a shallow 
cut-off, as described in §6.3 (Eqn. 2 with 
A brA; =670A, / = 2.8 and 5 = -1.6). The 
column is iV 2 o=l-0 and the dust mixture 
/di=0.6. 

Fig. 18. — A repetition of the spectrum of 
Pks 0232-04 of Fig. 9. In all models, the 
assumed near-UV index is a u =— 0.2, as in 
Fig. 9. The orange line corresponds to the 
standard intergalactic model introduced in 
§5, but combined with an intrinsic dust 
column of A^2o=0.18 with mixture /di=1-0. 
The cyan line (underneath the orange line) 
is the intrinsic dust model, but assuming 
that the powerlaw is modified by a shallow 
cut-off, as described in §6.3 (Eqn. 2 with 
A brfc =670A, / = 2.8 and 5 = -1.6). The 
column is A^" 2 o=0.62 and the dust mixture 
f D i=1.0 The blue line spectrum at the bot- 
tom corresponds to the class (D) spectrum 
of HE 1122-1649 (M /14 =0.44). The gray 
and magenta lines correspond to a power- 
law with a v = —0.6, multiplied by the cut- 
off function C\ of Eqn. 2, assuming a form 
factor / of 2.8 and 10, respectively. In both 
cases, the break wavelength is Afe r fc=670 A 
and the index change 5 = —2.0. 
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Fig. 19. — Black thin line: composite en- 
ergy distribution. Black dotted line: an 
SED consisting of a powerlaw of index 
a u —— 0.6 multiplied by the cut-off function 
C\ described by Eqn. 2. The parameters 
are the same as those inferred from the 
z = 2.8 transmission curve in Fig. 16, that 
is A ferfe =670A, 6=-1.6 and f=2.8. The 
darker gray line is a composite simulation, 
assuming intrinsic dust and the above pow- 
erlaw modified by the function C\. The 
dust column is A^ 2 o=0.8 with a mixture 
/di=0.3. Silver dashed line: the compos- 
ite simulation of previous Fig. 15, assum- 
ing intrinsic dust, but without an intrinsic 
cut-off (pure powerlaw). 
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